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Summary 

Adding material enhancements to a terahertz traveling-wave tube amplifier is investigated. Isotropic 
dielectrics, negative-index metamaterials, and anisotropic crystals are simulated, and plans to increase the 
efficiency of the device are discussed. Early results indicate that adding dielectric to the curved sections 
of the serpentine-shaped slow-wave circuit produce optimal changes in the cold-test characteristics of the 
device and a minimal drop in operating frequency. Additional results suggest that materials with 
simultaneously small relative permittivities and electrical conductivities are best suited for increasing the 
efficiency of the device. More research is required on the subject, and recommendations are given to 
determine the direction. 


Introduction 

Dielectrics have long been used to enhance the efficiency of waveguides and traveling-wave tube 
(TWT) amplifiers. The first field theory analysis of a dielectric-loaded helical TWT showed that greater 
electron beam coupling and higher efficiency were possible if a thin dielectric sheath was added to the 
inside of the beam tunnel (ref. 1). Additional studies found that an optimal dielectric constant existed for 
the beam tunnel sheath to better slow the primary mode and further increase electron beam interaction 
(ref. 2). Further, experiments have been done using ceramic enhancements to eliminate spurious modes in 
a gyrotron-traveling-wave tube amplifier and improve overall performance (ref. 3). In this paper, we will 
investigate applying various material enhancements to the waveguide surfaces of a proposed 0.4-THz 
traveling-wave tube amplifier to reduce circuit attenuation and increase efficiency. In addition to standard 
dielectrics, we will examine materials believed to display a negative refractive index (ref. 4), including 
anisotropic photonic crystals (ref. 5) and metamaterials with negative permittivities and permeabilities 
(refs. 6 to 9). 

Figure 1 displays a schematic of the amplifier under investigation. Pictured in blue is the slow-wave 
circuit, consisting of a serpentine-shaped folded waveguide with a center-mounted, square beam tunnel. It 
is modeled as a perfect vacuum. The yellow represents the walls of the slow-wave circuit, which are 
modeled as perfect electric conductors (PEC). In addition, a Cartesian coordinate system is selected such 
that the y-axis points down the center of the beam tunnel, the x-axis points horizontally across the width 
of the circuit, and the z-axis points vertically across the height of the circuit. Hereafter, we will refer to the 
y-, x-, and z-directions as the axial, horizontal, and vertical direction respectively. 

Models of the TWT were completed using COMSOL Multiphysics (ref. 10). Analysis was done using 
the RF Module's 3-D frequency-domain eigenmode solver. In order to efficiently simulate various 
dielectric parameters, the program's impedance boundary condition was utilized — a boundary effect that 
models the interaction between a material surface and a volume. In order to properly use this boundary 
condition, we make the assumption that the electric field will not penetrate deeply into the dielectric we 


NASA/TM— 2008-2 1 5059 


1 




Figure 1 . — Serpentine-shaped slow-wave circuit for terahertz 
traveling-wave tube amplifier. The model has yellow as perfect 
electric conductor and blue as vacuum. 

are modeling. This assumption is valid, considering the small thickness of the dielectric in relation to the 
width of the folded waveguide. 

The power output of a traveling-wave tube amplifier is strongly correlated to its cold-test 
characteristics, namely on-axis impedance and attenuation. On-axis impedance is a measure of the 
interaction between the electron beam and the electromagnetic wave inside the circuit. We calculate it for 
the /ith space harmonic by 


K„ 



(i) 


where Prf represents the RF power flow through the circuit (calculated from the axial Poynting vector), 
the /7th axial phase constant, and E n the magnitude of the on-axis electric field for the /7th harmonic 
(ref. 1 1). Likewise, attenuation is a measure of the resistive heating in a given circuit — that is, the power 
lost per cavity to the conducting walls. Given in decibels per cavity, we calculate circuit attenuation by 

a = 8.686-^- (2) 

2 Prf 

where L is the length of a single cavity and P L is the total power lost per unit length to resistive heating 
(ref. 12). We assume a conductivity for copper of a = 5.7xl0 7 S/m to calculate power loss. 

In order to reduce the number of cavities in a TWT amplifier and increase device efficiency, we must 
maximize the on-axis impedance and minimize attenuation. In order to more easily accomplish this, a 
single quantity was derived from curve-fitting simulation results with the NASA 2-Dimensional Coupled- 
Cavity Traveling- Wave Tube Code (ref. 13). It was found that the efficiency for our 0.4 THz TWT slow- 
wave circuit is proportional to the quantity: 

exp(0.3231 + 1.251- ^lK n - 1.532Va) (3) 

Methods to improve the performance of our device will center on optimizing this quantity. 

Our approach to increasing efficiency involves changing the boundary conditions inside the folded 
waveguide by simulating the presence of dielectric on its walls. Adding dielectric blocks the interaction 
between the electromagnetic wave and the conducting waveguide walls, thus reducing resistive losses and 
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decreasing attenuation. Furthermore, the orientation of the internal electric and magnetic field vectors 
change in the vicinity of the dielectric, disrupting the frequency of the simulated standing wave. 
Decreases in operating frequency, narrowing of bandwidth, and changes to on-axis impedance and 
efficiency are expected to occur. Our aim is to determine an optimal set of material parameters such that 
efficiency is maximized and operating frequency is kept high. In addition, we wish to determine a 
relationship between the cold-test characteristics of our TWT and the relative permeability and 
permittivity of our simulated dielectrics. 


Materials 

Three materials were investigated to increase the efficiency of our device: (1) the common dielectric 
silicon, with relative permittivity s r = 12.1, relative permeability p r = 1 , and electrical conductivity 
a = 10“ 12 S/m; (2) a general isotropic material with variable permittivity and permeability; and (3) a 
general anisotropic material with a variable permittivity tensor and p r = 1 . First, silicon was singled out 
because of its low conductivity and wide availability. It has found uses in vacuum electronics, 
waveguides, and traveling-wave tubes for decades, and its capabilities are well-understood. For those 
reasons, silicon represents an excellent starting point for our investigations. 

The second dielectric examined is a general isotropic material. This category includes most normal 
dielectrics — where the relative permittivity is positive and the relative permeability is unity — as well as 
negative-index metamaterials (NIMs) — where both permittivity and permeability are simultaneously 
negative. It should be pointed out that findings involving materials with negative material characteristics 
and non-unity permeabilities are largely theoretical results. In the case of metamaterials, which oftentimes 
consist of patterned or repeating geometries of both dielectrics and conductors, their macroscopic material 
parameters are not well-defined. In addition, the effect they have on incident electromagnetic waves is 
frequency-dependent. No such metamaterial to date has been constructed to operate at 0.4 THz, and few 
function over frequency bands wider than 2 to 3 GHz. Thus, present technology is poorly suited for our 
hypothetical amplifier. However, research is ongoing in the field of metamaterials. Dielectrics with 
tunable permittivities and permeabilities have been proposed, and work is underway to extend 
metamaterial properties to higher gigahertz frequencies. Consequently, as metamaterial technology 
matures, it may prove invaluable in the development of higher-power amplifiers and other vacuum 
electronics. 

The final material under investigation is a general anisotropic crystal. Commonly referred to as 
photonic crystals, they are periodic, latticed structures that have direction-dependent dielectric constants. 
Thus, an incident electromagnetic wave will experience different material properties depending on the 
direction and angle from which it enters the crystal. This property can be expressed mathematically by 
replacing each material constant with a rank 2 (two array indices are required to describe quantities) 
material tensor. For the sake of simplicity, we only consider anisotropies in the relative permittivity, with 
relative permeability set to unity. Moreover, we assume crystals that are oriented along a rectangular 
coordinate system. This zeroes out all permittivity values except those in the x-, y-, and z-directions, 
allowing us to express our tensor as 


8 = 


0 


"yy 

0 


0 

0 

^zz J 


( 4 ) 


By varying each of the permittivity components, a direction-dependent relationship can be found between 
the cold-test characteristics of our TWT and the material properties of our anisotropic dielectric. 

It is important to note that electrical conductivity has been removed from our discussion on isotropic 
and anisotropic materials. Instead, we assume a conductivity of silicon (a = 1CT 12 S/m) for all simulated 
material. We do this to ensure negligible dielectric losses from our material enhancements as well as to 
isolate the effects of permittivity and permeability on the cold-test characteristics. In the case of photonic 
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crystals, this could be argued as a valid assumption, given their uses in low-loss filtering and waveguiding 
(ref. 5). However, it becomes problematic when we apply this assumption to metamaterials. As mentioned 
above, their macroscopic material properties are not easily defined. In addition, since most metamaterials 
are composed of a dielectric base containing a repeating conductive geometry, they suffer from both 
resistive and absorptive losses (ref. 8). Further, Dewar calculated an effective conductivity of c e ff = 773 
S/m for a metamaterial consisting of a ferrimagnetic dielectric with an array of conducting wires (ref. 9). 
Therefore, assuming a low conductivity is not valid for current metamaterials. We retain the assumption 
for these simulations only to isolate the effects of permittivity and permeability. 

Methodology and Results 

Initial investigations involved using COMSOL's impedance boundary condition to simulate silicon 
(s r = 12.1, p r = 1, a = 10- 12 S/m) on various waveguide walls. Pictures were taken of the surface currents 
on each of the conducting boundaries and compared to the base model (fig. 2). In addition, dielectric 
losses, conductive losses, and attenuation were calculated to quantify the effects of adding silicon to 
various surfaces. Overall, our aim was to minimize attenuation by determining which walls would benefit 
most from the addition of dielectric. In the end, we found that placing silicon on the curved surfaces of 
the waveguide produced the largest drop in attenuation. 

These results are consistent with previous studies, which found conductive losses to be greatest along 
the inner curves of the serpentine waveguide (ref. 14). The reason for this is that the orientation of the 
internal magnetic field vectors — which lie tangent to the waveguide walls — change most rapidly in the 
vicinity of the inner curves. This rapid change produces substantial current densities, leading to greater 
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Figure 2. — Half of the slow-wave circuit is modeled. Contours represent 
surface current densities on the conducting waveguide walls. Note that 
density is largest along the inside curves. 
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Figure 3. — Half of the slow-wave circuit is modeled, with dark gray areas 
representing silicon boundaries. Contours are surface current densities. Largest 
densities appear near the inner curves of the waveguide at the interface between 
the silicon and the conductors. Magnitudes are larger than in figure 2. 

conductive losses in the circuit. Unfortunately, adding dielectric to the waveguide curves is not a perfect 
solution. The presence of a non-conducting material disrupts the electric and magnetic fields near the 
waveguide curves, causing surface currents to build up and intensify at the interface between the 
dielectric and the conducting waveguide walls (fig. 3). Ohmic losses are reduced on average with the 
addition of silicon, but this is mainly because fewer conducting walls are present to leach power from the 
internal wave. 

In addition to silicon, we tested to see how negative-index metamaterials (NIMs) would affect the 
cold-test characteristics of our TWT. Three simulations were conducted to test both silicon and an 
arbitrary NIM with relative permittivity s r = -1, relative permeability \i r = -1, and conductivity c = 10“ 12 
S/m. Two of the runs involved placing the materials on the inner curves of the serpentine waveguide, 
while the other involved only silicon on the outer curves. Eigenfrequencies, impedances, and attenuations 
were calculated for each, in addition to the relative efficiency calculated from (3). Figure 4 compares the 
values of each simulation to the cold-test characteristics of the non-enhanced TWT. 

As expected, placing dielectric on the inner curves of the waveguide produced the largest drop in 
attenuation at 0.4 THz. Unfortunately, it also led to huge losses in overall efficiency. For instance, silicon 
led to an attenuation drop of 45 percent from the original model, with a 60 percent drop in on-axis 
impedance and more than a 50 percent decrease in efficiency. Likewise, the NIM produced more than a 
55 percent drop in attenuation, at the expense of both a 75 percent loss in on-axis impedance and a 
65 percent drop in relative efficiency. On the other hand, enhancing the outside curves with silicon did 
not create as drastic a change. Indeed, the largest drop in operating frequency was only 2 percent across 
the entire dispersion curve. As a result, attenuation went down by 22 percent, on-axis impedance by 
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Phase shift per cavity, deg (b) Frequency, THz 



^ Frequency, THz 

Figure 4. — Cold-test characteristics for the base model compared 
to those with silicon on the inside waveguide curves, silicon on 
the outside curves, and negative-index metamaterial on the 
inside curves. Plots show (a) the dispersion curve, (b) the on- 
axis impedance, and (c) the attenuation. 


1 1 percent, and relative efficiency by 8 percent. Although these dielectric enhancements do not improve 
the cold-test characteristics of our TWT with its current dimensions, the results demonstrate that an 
appreciable drop in attenuation is possible — along with minimal changes to the efficiency of the device — 
if dielectric is placed only on the outside curves. Moreover, because the drop in operating frequency is 
small, only a small decrease in the circuit's width will be needed to shift the dispersion curve to a higher 
operating frequency. Figure 5 illustrates what the slow- wave circuit would look like with material 
enhancements on the outside curves. As before, the yellow structures are perfect electric conductors and 
the blue are vacuum; the red is dielectric. 

If we consider the fact that the NIM simulation produced the largest changes to the cold-test 
characteristics, it is conceivable that an optimal relative permittivity and permeability could exist to 
maximize the efficiency of our device. To test this hypothesis, parameter sweeps were conducted to 
determine the dependence of the cold-test characteristics on the relative permittivity and relative 
permeability. In general, it was found that, for a fixed relative permeability and a small electrical 
conductivity — with permittivity and permeability having opposite signs — both on-axis impedance and 
attenuation were found to increase as the permittivity approached zero. In addition, the two cold-test 
characteristics were found to increase at such a rate that the overall efficiency improved for small 
permittivities. Unfortunately, this result is not particularly realistic. Indeed, for such a material to exist — 
that is, to have both a small relative permittivity and a small conductivity simultaneously — it would have 
to absorb very little electric energy like a conductor yet resist the buildup of surface currents like a 
dielectric. No material in existence is known to display such hybridized properties. However, research in 
the field of metamaterials may one day produce a periodic structure with these properties. 
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Figure 5. — Serpentine-shaped slow-wave circuit with dielectric 
enhancements. Yellow areas are modeled as perfect electric conductor, 
blue as vacuum, and red as dielectric. 


In addition to negative-permittivity materials, we experimented with the addition of anisotropic 
dielectrics. Because of the serpentine shape of our waveguide, most cavity surfaces lay tangent to the 
horizontal x-direction. As a result, we only have significant electric field components along the axial and 
vertical directions inside the slow-wave circuit and not in the horizontal direction. Since no significant 
field exists to interact with the horizontal permittivity component, the permittivity tensor of any 
anisotropic material reduces to only ay-component (8 V> ) and a z-component (s zz ). Thus, all models tested 
remove the permittivity x-component (s xx = 0) and use s yy ± s zz . (The case where z yy = s zz represents a 
regular isotropic material, which we have already examined.) In addition, we set the relative permeability 
p r to unity and the electrical conductivity a to 10“ 12 S/m — a valid assumption given the relatively low-loss 
behavior of many photonic crystals. Both z yy and s zz were varied in a parameter sweep over positive and 
negative values, with operating frequency, on-axis impedance, attenuation, and efficiency calculated for 
each sweep. Multiphysics optimization algorithms were also run to determine which permittivity tensor 
produces maximum efficiency. 

In general, our anisotropic results are consistent with findings from the isotropic material tests and are 
similar for each of the following three cases: (1) both e yy and s zz positive, (2) both e yy and s zz negative, and 
(3) ^ positive and s zz negative. For each case, efficiency was found to increase as each tensor component 
approached zero. Figure 6 illustrates this trend graphically: 6a for case (1), 6b for case (2), and 6c for case 
(3). Note also that the largest efficiencies tend to occur when both and s zz are simultaneously small and 
equal. Thus, the optimal permittivity tensor for these three cases is not anisotropic at all. Instead, it 
confirms the result we found in our isotropic tests — that a dielectric-conductor hybrid with small relative 
permittivity and small electrical conductivity works best. Finally, scatter plots of efficiency versus 
frequency (fig. 7) for each case indicate that larger increases in the cold-test characteristics correlate with 
larger drops in the operating frequency. Consequently, finding an optimal efficiency would be 
problematic: if we increased efficiency, we would lower our operating frequency, forcing us to rescale 
our geometry. However, because our geometry has already been optimized for 0.4 THz, rescaling would 
in turn lower our efficiency. Not only that, we have no guarantee that this new higher efficiency would be 
any larger of an improvement over the existing model. Therefore, optimizing efficiency at the cost of 
large drops in operating frequency is not recommended. 
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s_yy > 0 and s_zz > 0 


s_yy < 0 and s_zz < 0 


(a) 




s_yy > 0 and s_zz < 0 



(c) 

Figure 6. — Device efficiency versus the y-component of permittivity for various anisotropic enhancements. 
The plots are for cases where the y- and z-components are (a) simultaneously positive, (b) simultaneously 
negative, and (c) the y-component positive and the z-component negative. 


Efficiency vs Operating Frequency 



Frequency, THz 

Figure 7. — Device efficiency is plotted against operating 
frequency for three cases of dielectric with anisotropic 
permittivity. Note the strong correlation between lower 
operating frequency and higher efficiency for these 
three cases. 


The fourth case (s ^ < 0 and s zz > 0) produces slightly different results. The cold-test characteristics 
are not as frequency-dependent as the other three cases. As the plot in figure 8 illustrates, there is not as 
strong a correlation between lower operating frequency and higher overall efficiency. Additionally, 
parameter sweeps were run to find an optimal permittivity tensor. It was determined that a small negative 
value for E yy and a value of around 3 for s zz produce the largest increase in overall efficiency. Figure 9 
presents the cold-test characteristics of the enhanced TWT, with s yy set to -0.0025 and s zz set to 3.1. 
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Although the device operates with both a higher efficiency and at a higher frequency than any 
combination of permittivities in the other three anisotropic cases, the frequency is still low compared to 
the original model, and substantial rescaling is again required. In addition, the bandwidth has narrowed 
from 200 GHz to less than 15 GHz — a problem that cannot be fixed by simply scaling the circuit 
geometry. 


Efficiency vs Operating Frequency for s_yy < 0 and s_zz > 0 



Figure 8. — Device efficiency is plotted against operating frequency 
for the anisotropic permittivity case of a negative y-component 
and a positive z-component. No correlation exists between 
higher efficiency and lower operating frequency. 


Comparison of Dispersion Curves 




0.30 0.35 


(b) 


0.40 0.45 0.50 

Frequency, THz 


0.55 0.60 



Frequency, THz 

Figure 9. — Cold-test characteristics for an amplifier enhanced with anisotropic-permittivity dielectric. 
Results are compared to the non-enhanced model in (a) dispersion curve, (b) on-axis impedance, and 
(c) attenuation. 
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Conclusions 


The incorporation of isotropic, negative-index, and anisotropic dielectrics into designs for a 0.4 THz 
traveling- wave tube amplifier has been investigated, and an optimal placement location has been found. 
Assuming a low electrical conductivity, the materials have been found to reduce conductive losses in the 
circuit. In addition, they have demonstrated sizeable reshaping of the circuit's internal electric and 
magnetic field vectors, particularly in regions that display large surface current densities. Unfortunately, 
nearly every dielectric enhancement investigated produced substantial decreases in both operating 
frequency and bandwidth, indicating that the dimensions of the TWT will need to be rescaled for 
upcoming simulations. In addition, the materials that demonstrated the greatest efficacy in increasing 
circuit efficiency are also the least understood and are currently unfeasible to obtain, construct, and 
incorporate. 

With regards to all enhancements modeled, anisotropic crystals were found to work no better than 
isotropic dielectrics to increase the efficiency of our device. Fabricating an anisotropic crystal to have 
properly oriented material tensors is a difficult task, especially on small scales. Moreover, simulating 
anisotropic parameters can lead to the excitation of spurious eigenmodes in finite-element codes, creating 
additional uncertainty in the model solution. Further, in most of the cases examined, the optimal 
anisotropic parameters obtained can be approximated by isotropic materials to achieve the same results. 
Therefore, no additional research can be recommended with regards to anisotropic material 
enhancements. 

Future work in this project will involve testing the dependence of the cold-test characteristics on the 
electrical conductivity of the material enhancements. Assuming a low conductivity may be valid for low- 
loss anisotropic structures like photonic crystals, but it produces unrealistic results for materials with 
small or negative relative permittivities. Silicon was used as a reference for electrical conductivity 
because of its negligible dielectric losses. Future models, however, will need to include these losses in the 
calculation of circuit attenuation to get a more complete view of using materials to enhance a traveling- 
wave tube amplifier. 

Now that a methodology is in place to simulate and test material enhancements, as well as trends 
relating a material's permittivity and permeability to circuit efficiency, additional tests will benefit from 
taking a more realistic look at existing dielectrics and metamaterials. Overall, as more information 
becomes available on the macroscopic material characteristics of metamaterials as well as their useful 
applications in a number of electron devices, a more accurate assessment on the results of this project can 
be made. 
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